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We investigate the dynamics of linear vortex lattices in anisotropic traps in two-dimensions and 
show that the interplay between the rotation and the anisotropy leads to a rich but highly regular 
dynamics. 
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INTRODUCTION 

Since their first observation in solid state supercon- 
ductors py] and liquid ^He [2[ arrays of quantized vor- 
tices have been one of the key features of systems ex- 
hibiting frictionless flow. More recently the observation 
of vortex lattices in gaseous Bose-Einstein condensates 
[3|, y has led to an re-newed interest in these systems, 
as ultracold atomic gases are highly configurable systems 
and allow for investigation of the lattices under a wide 
range of external parameters. Since then a large number 
of theoretical and experimental studies have been un- 
dertaken, which, for example, have examined the equi- 
librium properties of vortex lattices in different external 
potentials [5|, l6| or the dynamical properties of vortices 
lattices when disturbed [7|. At the same time new ideas 
for applications in quantum information [sl, l9| and preci- 
sion measurement [10| have been developed. 

The geometry of a vortex lattice is determined by the 
condition to minimize the energy of the overall system 
under the provision of fixed angular momentum. It there- 
fore depends on the interaction energy between the ve- 
locity fields of the vortices, which is a logarithmic func- 
tion of the distance between, and the underlying ex- 
ternal potential. For Bose-Einstein condensates carry- 
ing large angular momentum in isotropic harmonic traps 
this leads to the formation of a triangular Abrikosov lat- 
tices of vortices with winding number one y, lll|. For 
trap geometries that are not simply connected, such 
as a toroidal shapes, the system can also responds by 
forming persistent currents with higher winding num- 
bers [l2[. Other trap geometries and parameters, such as 



anisotropic traps [13 



or optical lattices [5| , no longer 
support an Abrikosov lattice, but instead form zigzag or 
linear, and square lattices respectively. 

While vortices and vortex lattices are interesting to 
study on their own merits 16|, |l7| , it has been suggested 
that they can be used in quantum information applica- 
tions through encoding of information into the winding 
number [8|, l9|. These numbers are so-called geometri- 
cal quantum numbers, which have the advantage of not 
suffering from any energetic instabilities and which are 
stable as long as the vortex cores are part of the system. 
Transferring the atoms into superposition states of these 
uantum numbers can be done either by optical methods 
or by making use of anisotropic external potentials 



|lq . However, due to the very dense energy spectrum 
of atomic superfluid systems it is a very difficult task to 
separate out exactly two flow states [l9|, l20| . Since the 
vortex-vortex interaction is dependent on both magni- 
tude and direction of rotation, however, one can in prin- 
ciple imagine a setup for a superfluid quantum computer 
analogous to the linear ion-trap model of Cirac and Zoller 
[2l| . In order to facilitate control and addressing of the 
individual vortex qubits it therefore would be preferable 
to reduce the number of nearest neighbours from six in 
the Abrikosov lattice, to, say, two in a linear crystal. It 
is known by today that two systems with which this can 
be achieved are waveguides IJ, |22| and anisotropic traps 



Here we consider the latter systems at medium rota- 
tion and with a small numbers of vortices. For currently 
typical experimental non-linearities it was recently found 
that even relatively small anisotropics, A ~ 2, lead to the 
vortices forming a linear crystals with equal spacing be- 



tween them [15|, |23|. However, the inhomogeneity of the 
background gas leads to a finite velocity of the vortex lat- 
tice within the cloud [2^ , which in the limit of large vor- 
ticity cancels out the external rotation [J]. For systems 
of medium vorticity in an anisotropic trap it is obvious 
that the linear shape of the crystal cannot be maintained 
during these dynamics and we will show here that this 
results in excitations within the background cloud which 
can be described by a standard hydrodynamics approach 
2 51 126| . More surprisingly, we also find that the excita- 
tions are oscillatory and do not destroy the ordering of 
the vortex lattice, even in the situation when the external 
rotation is switched off. 



In the following we will ffist briefly introduce our model 
and then examine the dynamics of the linear vortex lat- 
tices as a function of anisotropy by integrating the under- 
lying Gross-Pitaevskii equation for time scales of several 
hundred milliseconds. We determine the characteristic 
excitation frequencies of the vortex lattice and compare 
them to the known frequencies for anisotropic conden- 
sates. Finally, we remark on the stability of the vortices 
themselves. 



ANISOTROPIC TRAPPING POTENTIALS 



To simulate the dynamics of the vortices we consider a 
two-dimensional model of Bose condensed gas of particles 
of mass m in a harmonic trap, Vro = ^m(a;^x^ +a;^?/^), 
where the anisotropy is determined by the aspect ratio, 
A = ujy/ujx' The trap is rotating around the 2:-axis and 
we assume that a strong potential suppresses all relevant 
dynamics in this direction. For small, but non-negligible 
anisotropics (A ^ 2 for typical experimental systems), the 
ground state for medium rotation frequencies is given by 
a single line of vortices along the weak axis of the trap 
[15|. In the case of an odd number of vortices, one vortex 
will be located at the centre of the trap and the others 
will be symmetrically arranged around it. For an even 
number the same symmetry exists, however without the 
central vortex. Furthermore the vortex spacing remains 
roughly even along the soft axis due to the slow variation 
of the condensate density. 

For large rotational frequencies the condensate's efforts 
to distribute the vorticity as uniformly as possible lead to 
a self-consistent solid body rotation which cancels out the 
overall rotations of the lattice in the co-rotating frame. 
|24| . In the region of medium vorticity, however, in which 
only a few vortices are present, the concept of a diffused 
vorticity cannot be applied and the vortex lattice will 
try to carry out rotation in the co- moving frame. It is 
not surprising that in anisotropic traps this will lead to 
a dynamics that will have a feedback onto the density 
distribution of the atomic cloud. 

In the following we will consider these dynamics as a 
function of the trap anisotropy for the case of ^^Rb and 
currently realistic trapping frequencies. All our simu- 
lations are done by numerically integrating the Gross- 
Pitaevskii equation using a FFT/split-operator method. 



DYNAMICS AS FUNCTION OF ANISOTROPY 

The Hamiltonian for a Bose-Einstein condensate in an 
anisotropic trap and under rotation is given by 



H 



2m 



^HO - ^L, 



fl*p 



(1) 



where N is the number of atoms and Q is the frequency 
with which the trap rotates around the z-axis. As usual, 
for stability reasons Q is restricted to values smaller than 
the trapping frequency of the soft axis of the condensate. 
The two dimensional interaction strength is given by ^ = 



{4:7r?i^asc/m){muJz/27rhy^'^, with 



4.67 nm being 



the three dimensional scattering length for ^^Rb in the 
|2,1) state. 

We are interested in studying the behaviour of a 
ground state lattice as a function of time. While the 
behaviour of the condensate as a whole can be studied 




FIG. 1. The ground state density for a condensate of A/" = 10^ 
Rb atoms in a trap with A = 2.5 and cj^c = 10 x 27r Hz [upper 
panel] and the oscillatory trajectories of the vortices along the 
X-axis as a function of time [lower panel] . The centre vortex 
can be seen to remain stationary, while the other four vortices 
periodically oscillate away from their initial position. The 
dashed lines in the density plot indicate the range of move- 
ment in the ^-direction for each of the vortices (see Fig. [2|). 



using a h ydr odynamics approach in the Thomas-Fermi 
limit [251, l26|, the dynamics of the vortices themselves 
require an approach that is able to resolve the dynam- 
ics on the length scale of the healing length, such as the 
Gross-Pitaevskii equation. To have a clean initial state 
we therefore first find the ground state of the Hamilto- 
nian ([1]) by a numerical relaxation method and then in- 
tegrate this state in time while maintaining the rotation. 

The results of the numerical integration for a crystal of 
five vortices in a trap of aspect ratio A = 2.5 are shown in 
Fig. [TJ The upper half shows the starting configuration 
and the lower half displays the positions of the individual 
vortices as a function of time along the soft direction of 
the trap. The detailed displacement in the tight direction 
of the trap is shown in Fig. [2] and its range is indicated 
by the dashed lines in Fig. [TJ 

There are four immediate observations one can make 
from looking at these results: (i) the trajectories in both 
directions are symmetric for the pairs of vortices that are 
at equal distance from the centre of the trap. The cen- 
tral vortex remains stationary at all times; (ii) while the 
oscillation amplitudes for the two inner and the two out 
vortices are different, they are in phase with each other 
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FIG. 2. The position of the vortices in the ^/-direction corre- 
sponding to the situation shown in Fig. (TJ The top row shows 
the leftmost vortex and the bottom row the rightmost one. 



during the whole evolution; (iii) the vortex motion is pe- 
riodic and in particular the direction of rotation of the 
lattice reverses regularly and (iv) the paths of the vor- 
tices never cross and the linear ordering is maintained at 
all times. Closer examination shows that the trajectories 
of the two vortices on the right hand side and the two 
vortices on the left hand side follow the exact same pat- 
tern and are only different by a constant factor in their 
amplitude. This applies to the motion in the x as well as 
in the y direction and shows that the distance between 
neighbouring pairs of vortices is the same at any point 
in time. We have calculated these dynamics for different 
aspect ratios and found the behaviour consistently. As 
an example Fig. [3] shows the results for a trap with an as- 
pect ratio of A = 3.5, for which the same behaviour, apart 
from a smaller amplitude of the oscillations, is found. 

To better understand this pattern, we show in Table 
[J the dominating oscillation frequencies for the vortex 
oscillations in three different traps. For all anisotropics 
we find two dominating frequencies, uji and CJ2, and in 
the soft direction the lower frequency has a larger ampli- 
tude whereas in the tight direction the opposite is true. 
For increasing anisotropy the lower frequency decreases, 
whereas the higher frequency increases. This can clearly 
not be explained by straightforward oscillations in the 
harmonic potential, but indicates that the non-linearity 
of the sample plays a significant role. At the same time, it 
is also clear that the oscillation does not originate from 



FIG. 3. The ground state density for A = 3.5 [upper panel], 
and the position of the vortices along the x-axis [lower panel] . 
All system parameters are the same as in the A = 2.5 case 
and all observations made apply here as well with reduced 
amplitude. In the y-direction the movement is on the order 
of 1 X 10~^ m (not shown). 



a simply stationary vortex crystal inside an oscillating 
background cloud, as the cloud would be stationary in 
the co-rotating frame. However, we will show in the fol- 
lowing that both factors are playing an important role. 

Let us describe the dynamics of the background cloud 
using a hydrodynamic formalism and assume that the 
excitation of it are triggered and driven by the motion of 
the vortices. In fact, we assume that the rotation of the 
vortices leads to an infinitesimal rotation of the whole 
cloud and solve the classical hydrodynamical equations 
from there. This scenario was recently investigated by 
Lobo and Castin in order to find the scissors mode of a 
rotating condensate with vortices [26|. 

The density profile, p(r), and the velocity field, v(r), 
of a condensate in a rotating anisotropic trap in the 
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TABLE I. Frequencies dominating the oscillation of the vortex 
lattice from the numerical integration of the Gross-Pitaevskii 
equation (cji and 002) and oscillation frequencies from the hy- 
drodynamic model for the background cloud {ojhi and 0Jh2) 
for different values of the anisotropy. For comparison the 
trapping frequencies are shown as well. 



Thomas-Fermi limit in the co-rotating frame can be de- 
scribed by the equations of classical hydrodynamics as 

(St + V • V)v = V([/ -^pg)-2ftxv , 



m 

ft X {ft X r) 



dtp^V ' pv =0 



(2) 
(3) 



Because we have chosen the ground state with a station- 
ary linear vortex lattice as the initial condition for the 
cloud, the condensate is characterised by v ?^ for the 
velocity field and by the Thomas-Fermi solution for the 
density 

We assume that the rotation of the vortex lattice, and 
therefore the change in the associated flow, leads to a 
small, initial disturbance that corresponds to a rotation 
of the whole system with respect to the trap. It can 
be described by a small angle, 6>, from which the initial 
conditions for the linear perturbations to the above equa- 
tions are given as Jv = and 5p = Omxyiuo'^ — uOy)/ g. 
[26|. The subsequent dynamical evolution can then be 
calculated using the time-dependent polynomial ansatz 



5v = 5B{t)r 



r • SA{t)r] 



(5) 
(6) 



where SA{t) is a symmetric and 5B(t) is a general 2x2 
matrix and the chemical potential enters via the evolu- 
tion of the number c{t) [26| . The numerical solution of 
these equations shows an oscillation pattern which again 
is dominated by two main frequencies, uohi and U0h2-, and 
which we have tabulated in Table HI 

While these are not, and were not expected to be, a 
perfect match for the lattice oscillation frequencies, they 
are closer in value and in particular show the same trend 
of increase and decrease. This shows that the oscilla- 
tions of the lattice excite and drive oscillations within 
the background cloud, which in turn feeds back onto the 
lattice dynamics to create the continuous and periodic 
behaviour. This is also demonstrated in the fact that 
that the amplitude of the lattice oscillations is signifi- 
cantly larger than the amplitudes of the cloud. In Fig. |4] 
we show that the oscillation angle of the lattice cloud 
is clearly much larger than the one for the background 
cloud, which in fact stays very close to its initial dis- 
tribution. Nevertheless, the whole dynamics is highly 
predictable and in particular the order of the vortices is 
maintained throughout. 




FIG. 4. (Colour online) Density of the cloud at t = 0.0072s 
[upper panel] and t = 0.0108s [lower panel]. The full white 
line indicates the tilted axis of the vortex crystal and the bro- 
ken white line the horizontal starting position. The oscillation 
of amplitude of the background cloud is too small to be visible 
in this plot. 



in the following investigate the dynamics of the cloud 
in the absence of external rotation. From the results 
above it should be clear that the background flow will be 
strongly disturbed by the now non-rotating trapping po- 
tential and one might expect the linear vortex crystal to 
break up into a random pattern. As the system is finite, 
it is clear, though, that after some evolution the linear 
form has to be re-assumed. 

In Fig. [5] we show an example of the time evolution 
of a cloud with initially five vortices aligned along the 
soft axis (a) and an aspect ratio of A = 2.5. The vor- 
tex crystal can be seen to rotate against the cloud, still 
maintaining its linear shape after t = 0.03 s (b). After 
a further 0.03 s the linear geometry is slightly disturbed 
when the vortices are passing through the narrow direc- 
tion of the condensate (c), but they finally assume the 
linear shape again after about 0.1 s (d). As for the case 
of maintained trap rotation, the outer most pair of vor- 
tices follows a symmetrical path, as do the inner pair 
of vortices. However once the narrow section is passed 
the pairs on each side have changed their order for the 
first time(see Fig. [6]) and this process will be repeated 
over and over again during the ensuing evolution. The 
vortices also regularly form a linear crystal again (see 
Fig. [Sl^d)), however not necessarily along the horizon- 
tal axis. This behaviour persists for higher ratios of the 
anisotropy, until the condensate itself decays too fast, 
due to the anisotropic barrier it runs into, to be able to 
define a clear vortex crystal. 



DYNAMICS IN THE ABSENCE OF ROTATION 

It is interesting to investigate the robustness of the 
order of the vortex lattice under even stronger pertur- 
bations of the background cloud and we will therefore 



STABILITY OF VORTICES 

Let us finally briefly comment on the stability of vor- 
tices in the lattice. As the systems we investigate are 




x(|jm) 



FIG. 5. (Colour online) Dynamical evolution of the density 
and phase for a trap of aspect ration A = 2.5 for which the 
rotation has been switched off. Panel (a) shows the initial 
linear crystal at t = s and the subsequent panels show the 
times t = 0.03 s (b), t = 0.06 s (c) and t = 0.1 s (d). Panel (c) 
shows the disturbed crystal as outer vortices pass the narrow 
direction of the condensate as they rotate around the cen- 
tral vortex and panel (d) shows the reformation of the linear 
crystal. At this time the left and right pairs of vortices have 
rotated around the central vortex, and the inner and outer 
vortices on each side have swapped position. The phase plots 
show that the charge of the vortices is maintained throughout 
the time evolution. 



not rotationally symmetric, they do not need to conserve 
angular momentum and it is possible for the vortices to 
reverse their direction of flow [27|,|28|. However, this tun- 
neling between different rotational states is suppressed at 
larger interaction strengths and as all our simulations are 
carried out in the Thomas- Fermi limit, we have checked 
that all vortices keep their initial rotational direction 
throughout the simulations. To demonstrate this for one 
example, we show the phase distributions for the case of 
the non-rotating cloud in Fig. [6l 



CONCLUSION 

We have numerically investigated the dynamics of a 
linear vortex lattice in an anisotropic, two-dimensional 
Bose-Einstein condensate. Starting from the ground 
state of the Gross-Pitaevskii equation, we have found 
that the vortices retain an ordered structure both in the 
presence and absence of rotation for time-scales on the 
order of 1 s. If the rotation of the cloud is maintained, 
the vortex lattice remains linear and only carries out well 
defined oscillations within the background cloud. With 
increased anisotropy oscillation amplitude decreases until 
finally the crystal does not move within the cloud. This 
allows for the use of the vortex lattices in applications 




FIG. 6. Position of the vortices along the x-direction for Q = 
and for A = 2.5. The values for the initially outer pair of vor- 
tices are given by the full line and the inner pair of vortices 
by the broken. It can be seen that the paths remain symmet- 
rical for the respective pairs, however they switch positions 
periodically. 



such as, for example, in quantum information. In the 
absence of rotation we found that the outer vortices cir- 
cle around the central vortex, periodically reversing their 
order and the lattice and linear crystals are periodically 
reformed for not too large anisotropics. Once these get 
too strong, the background cloud gets too distorted to 
still talk about a well defined vortex lattice. 

The work was supported by Science Foundation Ire- 
land under project number 05/IN/I852. 
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